Global Analysis of Synchronization in Coupled Maps 
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We introduce a new method for determining the global stability of synchronization in systems 
of coupled identical maps. The method is based on the study of invariant measures. Besides the 
simplest non-trivial example, namely two symmetrically coupled tent maps, we also treat the case 
of two asymmetrically coupled tent maps as well as a globally coupled network. Our main result 
is the identification of the precise value of the coupling parameter where the synchronizing and 
desynchronizing transitions take place. 

PACS numbers: 05.45.Ra, 05.45.Xt 



INTRODUCTION 

Although the phenomenon of synchronization of non- 
linear oscillators was discovered already in 1665 by 
Huyghens, it has received systematic attention only re- 
cently P, , owing mainly to newly found applications 
and the understanding of nonlinear systems that we 
have achieved through modern methods. The field re- 
ceived an impetus after it was discovered by Fujisaka- 
Yamada Q and Pikovsky [|| and further analyzed by 
Pecor a- Carroll Q that even chaotic trajectories can syn- 
chronize. Since then the phenomenon has found applica- 
tions in secure communication, coupled Josephson junc- 
tions arrays etc. From a completely different perspective, 
it was observed 0,0 experimentally that synchronous fir- 
ing of neurons was important in feature binding in neural 
networks. Distant groups of neurons rapidly synchro- 
nize and desynchronize as the input stimulus changes, 
each synchronized group representing a collection of fea- 
tures belonging to the same percept. Thus a systematic 
understanding of synchronization and desynchronization 
has become important in neural information process- 
ingiilS 

As a result of these developments various coupled dy- 
namical systems have been studied that are either con- 
tinuous or discrete in time, either continuously coupled 
or pulse coupled etc. and different types of synchroniza- 
tions have arisen such as phase synchronization or com- 
plete synchronization, etc. Coupled maps E3 have been 
used to a grea t extent to understand synchronization in 
particular |l2L 

EE D El E3, E3, EM Em EE B and 

complexity in dynamical systems \22L 1 2.1 124| in general. 
The paradigm here consists of identical individual maps, 
typically iterates of a functional equation like the logistic 
one that produce chaotic dynamics. These maps then 
are coupled, that is each of them computes its next state 
not only on the basis of its own present state but also 
on the basis of those other ones that it is coupled to. 
With appropriate conditions on the coupling strengths, 
the individual solution is also a solution of the collective 
dynamics, that is, when all the maps follow their own 
intrinsic dynamics synchronously, we also have a solu- 



tion of the coupled dynamics. This is the simplest case 
of synchronization. That synchronized collective dynam- 
ics, however, need not be stable even if the individual 
ones are. Therefore, the study of synchronization essen- 
tially depends on a stability analysis. One may use lin- 
ear stability analysis or global stability analysis (for some 
recent formulations, see HEEBIES) ^ n the linear sta- 
bility analysis one assumes that the synchronized state 
follows the same dynamics as the individual uncoupled 
dynamics and then studies the stability of the dynam- 
ics transverse to the synchronizing manifold. However 
it has been observed |2j| in coupled systems with delays 
that the synchronized dynamics can be quite different 
from the individual uncoupled dynamics. Also, it has 
been demonstrated |29l |30| that the linear stability anal- 
ysis can actually fail in a nonlinear setting, due to global 
nonlinear effects. For a global stability analysis, one has 
to guess a Lyapunov function which then gives a criterion 
for synchronization. Thus the linear stability takes only 
the local dynamics into account and it makes a restricting 
assumption on the synchronized dynamics. In contrast, 
the global stability analysis, while global, depends on the 
choice of the Lyapunov function and may not always give 
optimal criteria. In general optimal global results are dif- 
ficult to obtain. In many cases numerical results indicate 
the global stability of the synchronized state, but it is dif- 
ficult to prove stability rigorously. Thus it is necessary 
to develop a method to study synchronization which goes 
beyond these drawbacks. That is, it should be global 
and should not make any reference to the synchronized 
dynamics. This is what we achieve in this paper. We 
present a new global stability method that leads to sharp 
results, that is, it can determine the critical value of the 
coupling parameter above which global synchronization 
occurs. In order to exhibit the principal features, we first 
apply the method to the simplest non-trivial case, namely 
two symmetrically coupled tent maps. Then we consider 
asymmetric coupling and a globally coupled network. 

Our method uses invariant measures or stationary 
densities|3lj to study synchronization. A stationary den- 
sity is the density obtained by starting from some initial 
distribution of the initial conditions and letting the sys- 
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tern evolve for a long time. Though the complete char- 
acterization of such densities can be difficult, we shall 
demonstrate below that it may suffice to study its sup- 
port to understand synchronization and this can be eas- 
ier than studying the whole measure. Clearly, this ap- 
proach is inherently global. The key idea is to study 
the support of the invariant measure for different cou- 
pling strengths and to determine the critical value of the 
coupling strength when the support shrinks to the syn- 
chronized manifold. We consider the following coupled 
dynamical system 



An 



Af(X n ) =: S(X n ) 



(1) 



where X = (xi, X2, • ■ ■ , xn) t is an iV-dim. column vec- 
tor, A is an N x N coupling matrix and / is a map from 
il = [0, 1] N onto itself. As already mentioned, in the 
present paper we consider the tent map defined as: 



, . _ J 2x < x < 1/2 
I( - X) ~ \ 2-2x l/2<x<l. 



(2) 



Different coupling matrices are chosen. 



obtain a singular measure asymptotically. However, if, 
for example, the map is expanding everywhere then an 
SRB measure is absolutely continuous w.r.t. to Lebesgue 
measure. Therefore here we restrict ourselves to one par- 
ticular map that is expanding, namely the tent map. We 
want to investigate how such a density depends on the 
coupling strength and would like to see when its support 
shrinks to the synchronization manifold. 

There are various ways to find invariant measures. One 
approach is to use the so called Frobenius-Perron opera- 
tor. It is defined as: 



Pp{X')d N X' = / / p{X')dX' (4) 

J JS- 1 (D) 



If we choose D — [0, xi] x • • • x [0, xjy] then we get 

p(X')dX' (5) 



Our S is not invertible. In fact, it has 2 N disjoint parts. 
Let us denote them by S" 1 , i = 1, 2 N . If X £ Q, since 
/ is symmetric, we get 



INVARIANT MEASURES 

A measure p is said to be invariant under a transfor- 
mation S if p(S~ 1 (A)) — p(A) for any measurable subset 
A of f2. For example, a Dirac measure concentrated at 
a fixed point of the transformation is clearly invariant. 
An invariant measure need not be unique as can be seen 
immediately from the fact that if there are several fixed 
points then the different Dirac measures at these fixed 
points as well as their linear combinations are all invari- 
ant. A natural way to obtain invariant measures is to 
simply iterate any measure under the transformation S 
and take an asymptotic (weak-*) limit of these iterates. 
A sequence p n of measures converges in the weak-* sense 
to a measure p when for all continuous functions / on 



lim / f(x)p n (dx) 
In 



f(x) p{dx) 



(3) 



The usefulness of this concept derives from the fact that 
bounded sets are weak-* precompact, that is, contain a 
weak-* convergent subsequence, see e.g. j^. In particu- 
lar, this applies to the dynamical iterates of some initial 
measure. For example, one may start with any Dirac 
measure. Our interest, however, is not in such singular 
measures as they do not sample the whole phase space. 
We would like to start from a distribution of initial condi- 
tions spread over the whole phase space (say uniformly) 
and study its evolution and what kind of asymptotic limit 
it leads to. This is the idea of the SRB-measures, as they 
are called after Sinai, Ruelle, and Bowen. For some dy- 
namical systems, any invariant measure is singular. In 
such cases even if we start with a uniform density we 



Pp(X) = J-\X)^2p(S^(X)) 



(6) 



i=i 



where J" 1 (A) = \dS- 1 (X)/dX\ 



SYMMETRICALLY COUPLED MAPS 

In this section we choose a symmetric dissipative cou- 
pling given by 



A 



a — e e 
e a — e 



(7) 



where < e < 1 is the coupling strength and < a < 1 
is a parameter. This choice for A also satisfies the con- 
straint that the row sum is a constant. This guarantees 
that the synchronized solution exists. 



The case a — 1 

With this choice of A we get the following functional 
equation for the density. 



Pp(x,y)- 



1 



r\p(px/2-<ry/2,->Yx/2 + py/2) 



4|l-2e| 

+p(l - Px/2 + 7 y/2, -jx/2 + (3y/2) 

+p(f3x/2 - 7J//2, 1 + -/x/2 - Py/2) 

+p(l - Px/2 + jy/2, 1 + 1X /2 - Py/2)] (8) 

where 7 = e/1 — 2e and P = 1 + 7. Since we know that a 
point belonging to f2 does not leave f2, all the arguments 
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FIG. 1: The evolution of the support of the invariant measure 
obtained by starting from the Lebesgue measure supported on 
the square [0, 1] x [0, 1] with e < 1/4. a) The parallelogram 
ABCD obtained after the first application of A b) The rhom- 
bus AFIE obtained by applying / to the rectangle ABCD 
in a. c) The rhombus AFIE one gets when A is applied to 
the rhombus in b. d) The next application of / leads to the 
rhombus AMNL. 



of p on the right hand side of the above equation should 
be between and 1. This gives us four lines < fix/ 2 — 
7Z//2 < 1 and < —jx/2 + /3y/2 < 1 which bound an 
area, say T. The support of the invariant measure should 
be contained in T n £1. 

In order to study the evolution of supp p we study the 
iterates of Vt under S — Af. / maps f2 onto f2. Applica- 
tion of A then leads to the parallelogram shown in Fig.[TJi. 
The next application of / leads to Fig. ^p. We use the 
following convention for labelling the points throughout. 
A point is labelled by the same letter in its image. As a 
result, some points in Fig.^3 have two labels. We note 
that the boundary in Fig. is completely determined 
by the part of the area in the upper right quadrant in 
Fig. [IJi. This is essentially because l(GI) > 1{EI) and 
I (HI) > l(FI). Let us call this condition A. So the areas 
in other quadrants get mapped inside the image of this 
area. Now we operate A again to obtain Fig. ^ and / 
again to get Fig.^Jl. Here we have assumed that e is suffi- 
ciently small so that the x-coordinate of F in Fig.^ is less 
than 1/2 and by symmetry the y-coordinate of E is less 
than 1/2 (condition B). One can notice that, although the 
internal structure of the density becomes more complex, 
the boundary is the same as that in Fig. ^p. Also, by 
continuing this procedure one can see that the quadrilat- 
eral AMNL in Fig.^i is the support of the density which 
remains invariant provided e is sufficiently small |40j. 



1 -e e 
e l-e 



l-e 
1 



2e 

-2c+2c 2 



(9) 



As a result, the condition B leads to e < 1/4. So when e 
satisfies this condition we get the above area as an invari- 
ant area. In |3^|, Pikovsky and Grassberger have made 
a similar observation for an asymmetric generalisation of 
the tent map and later Glendenning |3^] has carried for- 
ward the analysis mathematically. However there is an 
important difference between the work in piij l and our 
work. The aim of the ref. j^H is to study the blowout 
bifurcation and the Milnor attractor. As a result, there 
the analysis is restricted to the range of coupling param- 
eters between the blowout bifurcation transition and the 
complete synchronization transition (this range is absent 
in the case of the symmetric tent map), whereas we are 
studying the transition to complete synchronization. We 
carry out the complete analysis using purely geometric 
arguments preserving the global nature of the result. In 
the following two sections we shall also extend the anal- 
ysis to more general situations. 

We remark that the procedure carried out so far is 
applicable to a general class of nonlinear maps. Thus, 
by this method, one immediately obtains a value of e 
below which there is no synchronization provided one can 
argue that there exists an absolutely continuous invariant 
measure. Further discussion of this point is postponed to 
the concluding section. 

Now we have to check what happens for e > 1/4. So 
we carry out the iterations with this condition in mind. 
As a result in place of Fig. QJ: we obtain Fig. 01. Now 
if we apply / we get Fig. |2h>. We see that the boundary 
AEMNLF we obtain now is slightly different from that 
in Fig.^i, the corners having been chipped off. The next 
application of A gives us Fig- Et- It can be easily checked 
explicitly that the point E lies on the left and the point 
M lies on the right of the line x = 1/2 for any e greater 
than 1/4 j^jj. This implies that the boundary obtained 
by the next application of / (Fig. 01) looks similar, that 
is, the number of vertices remains the same. In fact, 
except for the slanted portions it is exactly the same and 
we set out to argue that it is indeed different on those 
small portions. That is we cannot have a set with such a 
boundary as an invariant area. 

In order to do this, let us first note that the magnitude 
of the slope of any line of the form y = mx + c does not 
change with / as it scales both axes by the same factor. 
And A maps such a line to a new line given by 



y 



+ (1 — e)m ec(e + (1 — e)m) 



1 



1 — e + me 

Now if the slanted portion is stable 
shouldn't change under A, so we have 

e+(l 



+ (1- 
then 



-e)c (10) 
its slope 



l-e 



m 



(11) 




A A,N 



FIG. 2: The evolution of the support of the invariant measure 
obtained by starting from the Lebesgue measure supported on 
the square [0, 1] x [0, 1] with e > 1/4. a) The rhombus AFIE 
obtained by applying A to the rhombus AFIE in Fig. 
b) The polygon AEMNLF one gets after the application of 
/ to the rhombus AFIE of a. c) The next application of A 
transforms the polygon of b to this polygon AEMNLF. d) The 
new polygon ALOQPM obtained by the next application of 
/ to the polygon in c. 

This implies m 2 = 1 (m = 1 in our case since we are in- 
terested in the region where both coordinates are greater 
than 1/2 and / does not change the slope in this region). 
So the slanted portion in the Fig. if it is stable, has 
to be of the form y = x — (1 — 2e)c. / changes this in- 
tercept to —2(1 — 2e)c and for it to be stable this should 
be equal to — c owing to the symmetry. This implies that 
either 6 = 0ore = l/4. Therefore such a slanted por- 
tion cannot be stable for e > 1/4 and the whole area 
cannot be stable either. This shows that e c = 1/4 is the 
critical value below which there exists an invariant area 
and hence no synchronization. When e crosses this e c the 
area becomes unstable and collapses the synchronization 
manifold to the line x — y and we get synchronization. 

If the condition A is violated as well in Fig. during 
further iterations then in that case the boundary would 
look something like the one shown in Fig. instead of 
that in Fig. Hi. But since LMNL > LEAF, there is 
always a part of nonzero length of the segment MP in 
Fig. [St which appears as a boundary in Fig. (segment 
MT). And then the above argument will apply to this 
segment. 

By symmetry it can be argued that this area again 
becomes stable if we increase e further at e = 3/4 and 
this leads to desynchronization. Needless to say that the 
result agrees with that obtained by the linear stability 
analysis. 

It is interesting to note that the transitions of synchro- 



A.N 

FIG. 3: This polygon would be obtained as the support of 
the evolving measure during the course of iterations if the 
condition A (see the text) was violated. 

nization and desynchronization are discontinuous in the 
sense that the area of the support abruptly changes value 
at the critical thresholds. 

Finally, we end this section by making a couple of ob- 
servations which might be useful in future refinements 
and/or generalizations of this analysis to other maps. 
Firstly, it is interesting to note that the vertex F in 
Fig ^3 or c becomes a periodic point of period 2 ex- 
actly at e = 1/4 and it looses its stability for larger e. 
Secondly, for e < 1/4, the boundary of the invariant area 
is generated from the boundary of f2 and hence is sta- 
ble whereas for e > 1/4 there is a small portion of the 
boundary (segment EM in the Fig. [21 b or c) which does 
not satisfy this property and becomes unstable. 

The case of general a 

In this case the synchronized dynamics is different from 
that of the individual maps. This can be seen by taking 
the initial point as Xq = (z, z) which lies on the syn- 
chronization manifold which then leads to the dynamics 
governed by X n+1 = af(X n ). 

Now we start the iterations of fl. The first application 
of A leads to the Fig. 0Ji and then / maps it to Fig. 0Jd. 
One can see that if e is sufficiently small so that the 
point H remains on the left of the line x = 1/2 after 
application of A then the IHCG remains the same in the 
subsequent iterations whereas the LFAE decreases with 
the iterations. This leads to the Fig. [S] asymptotically. 

The above condition on e can be calculated as for the 
case a = 1 and we get e < (2a — l)/4. By following 
similar arguments as for the a — 1 case it can be shown 
that this area is unstable for e > (2a — 1)/4 and it reduces 
to the line x — y. Again by symmetry it can be seen 
that the desynchronizing transition will take place at e = 
(2a + l)/4. 
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FIG. 4: The evolution of the support of the measure for the 
coupling matrix A with general a. The iterates are started 
from the square [0, 1] x [0, 1] and e < (2a - l)/4. a) The 
parallelogram ABCD obtained after the first application of 
A. b) The next application of / transforms the parallelogram 
in a to this figure. 



FIG. 5: The invariant area CHIG obtained asymptotically 
for the evolution with the coupling matrix A with general a. 



As a result, £ (61,62) and 5'(e2,ei) synchronize simulta- 
neously when the values of e\ and 62 are varied. 

In fact, much more is true. Let us define an operator 



St 



S( 



£1 + te 2 tei + e 2 , 
l + t ' t+1 ' 



S(ei,e 2 ) + tS(e 2 ,e 1 
l + t 



■15) 



In particular Sq — S(ei,e2) and <Soo = 5(62,61). Now 
it can be seen that StX lies on a straight line between 
SqX and SooX and all the three points lie on the same 
side of the line x — y. This in particular implies that, for 
given values of ei and 62, the boundary of the support of 
the invariant measure of St lies between that of SqX and 

x 



SqqX. In fact, if we define D 



\x — y\ then we see 



immediately that D(St 



\l- ei -e 2 \\f(x)-f(y)\ 



is independent of t and that, consequently, as ei and e 2 
are varied St synchronizes independently of t. 

Now when t = 1 we have Si = S{^^-, which 
is symmetric. So we can use the result of the previous 
section which tells us that the system synchronizes when 
ei + 62 = 1/2 and desynchronizes when t\ + e 2 = 3/2. 

It should be noted that the invariant measures gener- 
ated by StS for different t are not the same. But since we 
are interested in the synchronization property, i.e., the 
value of e at which the support of the measure shrinks 
to the line x = y, we use the identity l|15|) which gives 
another equivalent symmetric dynamics. And since this 
equivalent dynamics is symmetric we know its synchro- 
nization properties. 



GLOBALLY COUPLED MAP NETWORK 



ASYMMETRICALLY COUPLED MAPS 

In this section we again choose N = 2 but now 
1 - ei ei 



A = 



£2 1 - £2 



(12) 



where < £1,62 < 1. We write S(ei,e2) = Af. We can 
immediately see that 



S {ll . <2 )[ X y \ =RS(e 2 ,e 1 )( y x 



(13) 



where 



R 



1 
10/' 



i.e., the system remains unchanged under the transforma- 
tion ei <-> e 2 and x <-> y. This implies that, if p ei>e2 {x,y) 
and p e2l€1 (x, y) are the stationary densities generated by 
5(61,62) and S(e 2 ,ei) respectively, 



Pe lt e a (x,y) = p e2<ei (y,x). 



(14) 



Now we extend our analysis to a higher dimensional 
case. We consider a globally coupled network of N maps 
with all couplings of equal strength. As a result 



.4 



( 1 £ 

e/(N-l) 



e/(N-l) 
1-e 



e/(N-l)\ 
e/(N-l) 



\e/(N-l) e/(N-l) 



1 



(16) 



/ 



Now we consider the dynamics in the 2-dim space 
spanned by ei and (1,1, ■•■,!)■ By symmetry this dy- 
namics is the same as that in any other similar subspace, 
i.e., spanned by and (1, 1, ■ ■ ■ ,1). And all of them 
synchronize simultaneously. We take the orthogonal vec- 
tors ei and (0, 1, • • • , 1) and a general vector in this space 
X = (x,y,y, ■ ■ ■ ,y) T . Now 



AX = 



( (l-e)x + ey ^ 

jv— 1 x (1 — N—i )y 



\ w=jx + (1 - 7 ^ I )y j 



(17) 
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As a result, the reduced dynamics of our system in this 
2-dim subspace is given by an effective coupling matrix 



.4 



eff 



l-e e 
T (1 



N-l' 



(18) 



and from the previous section we know that this dynam- 
ics synchronizes when e + e/(N — 1) = 1/2, i.e., e = 4^". 
This result again agrees with the one obtained by the lin- 
ear stability analysis [35l | . 



CONCLUDING DISCUSSION 

In this paper we have proposed a method that uses in- 
variant measures and in particular their support to study 
synchronization. We have demonstrated it on one exam- 
ple, namely the tent map where we were able to study 
the stability of the support of the invariant measure us- 
ing purely geometric arguments. In the future, different 
techniques to deal with invariant measures and their sup- 
port will have to be developed in order for this method 
to be of wide use. 

We should note that the notion of convergence we have 
been studying here, namely convergence of the support 
of the invariant measure as encoded in the convergence 
of its boundary in general is different from the weak-* 
convergence introduced above. That is, there exists the 
possibility that the measure can become zero asymptoti- 
cally everywhere except on the synchronization manifold 
without making the boundary unstable. In other words, 
the measure can become singular. Here, we have consid- 
ered the tent map which is expanding everywhere. But 
it is not sufficient to have the individual map expanding, 
rather it should be the coupled map that should guaran- 
tee the existence of the absolutely continuous measure. 
In fact, it turns out that, in all the cases considered here, 
the combined map ceases to be expanding exactly at the 
synchronization threshold. Clearly this is a consequence 
of the everywhere expanding property of the map chosen. 
While this guaranteed for us the existence of the abso- 
lutely continuous measure for e < e c and validated our 
use of the present method, it also raises a question as to 
what happens in the case of other maps. In fact, there 
exists a considerable activity in Ergodic Theory concern- 
ing the existence of an absolutely continuous invariant 
measures for various discrete dynamical systems (see for 
example [13, H3, So for general mappings our ap- 

proach will have to be complemented by the results in 
this field. 

Other than considering different maps one can also try 
to apply this method to more complicated connection 
topologies. Various networks, such as, random, scale- 
free, small world etc. can be considered. A characteri- 
zation of complete measures is also needed. We have al- 
ready taken a step in this direction |39j |. It is important 
to note that the knowledge of the complete measure was 



not necessary to study the synchronization since we used 
only the support of the measure. However, a method 
to find out that support is the minimum prerequisite for 
this approach to be useful. 

One of us (KMK) would like to thank the Alexander- 
von-Humboldt-Stiftung for financial support. 
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